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ABSTRACT 

Recently, the HESS collaboration has reported the detection of 7 -ray emission above a few hundred GeV 
from eight new sources located close to the Galactic Plane. The source HESS J1813-178 has sparked par¬ 
ticular interest, as subsequent radio observations imply an association with SNR G12.82-0.02. Triggered 
by the detection in VHE 7 -rays, a positionally coincident source has also been found in INTEGRAL and 
ASCA data. In this Letter we present MAGIC observations of HESS J1813-178, resulting in the detec¬ 
tion of a differential 7 -ray flux consistent with a hard-slope power law, described as dN 7 / (dAdtdE) = 

(3.3 ± 0.5) x 10 _ 12 (£'/TeV ) _21±0,2 cm _ 2 s _ 1 TeV -1 . We briefly discuss the observational technique used, 
the procedure implemented for the data analysis, and put this detection in the perspective of multifrequency 
observations. 

Subject headings: gamma rays: observations, supernova remnants: individual (HESS J 1817-178, G12.82- 
0 . 02 ), acceleration of particles 


1. INTRODUCTION 

In the Galactic Plane scan performed with the HESS 
Cherenkov array in 2004, with a flux sensitivity of 3% Crab 
units for Y-rays above 100 GeV , eight sources were discov¬ 
ered (lAharonian et ~af] 12005 aihl) . One of the newly detected 
7 -ray sources was HESS J1813-178. At the beginning, HESS 
J1813-178 could not be identified and was assumed to be a 
“dark particle accelerator”, without reported counterparts at 
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lower frequencies. 

Since the original discovery, HESS J1813-178 has been 
associated with the supernova remnant SNR G12.82-0.0 2 

dUbertini et alJ2005tlBrogan et alfctMlHelfand et al.l2005l) . 

One may still not exclude this_coincidence being the result 
of just a chance association. lAharonian et all ( 2005ah state a 
probability of ( 6 %) that one of their new sources is by chance 
spatially consistent with a SNR. Nevertheless, the properties 
of SNR G12.82-0.02, the multifrequency spectral energy dis¬ 
tribution, and the flux and spectrum of the high energy 7 -rays 
detected from this direction appear to be consistent with a 
SNR origin. 

HESS J1813-178 has been found to be n early point like 
(it source = 2.2') bv Aharonian et al J il2005bl) . Given the size 
of the SNR, the angular resolution of the HESS telescope, and 
the depth of the observations, the source size does not rule out 
a possible shell origin. The 7 -ray source lies at 10' from the 
center of the radio source W33. This patch of the sky is highly 
obscured and has indications of being a recent star formation 
region (IChurchwelli1990t) . 

In this Letter we present Major Atmospheric Gamma Imag¬ 
ing Cherenkov telescope (MAGIC) observations of HESS 
J1813-178. We briefly discuss the observational technique 
used and the procedure implemented for the data analysis. We 
derive a high energy 7 -ray spectrum of the source, and ana¬ 
lyze it in the perspective of multifrequency observations. 
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Fig. 1.— Sky map of 7 -ray candidate events (background subtracted) in 
the directions of HESS J1813-178 for SIZE > 600 ph. el. (corresponding to 
an energy threshold of about 1 TeV). Overlayed are c ontours of 90 cm VLA 
radio (black) and ASCA X-ray data (green) from IBrogan et alJ <20051) . The 
two white stars denote the tracking positions W1,W2 in the wobble mode. 

2. OBSERVATIONS 

MAGIC (see e.g., iBaixeras et alJ (120041) : ICortina et alJ 
< 2001 ) is currently the largest single dish Imaging Air 
Cherenkov Telescope (IACT) in operation. Located on the 
Canary Island La Palma (28.8°N, 17.8°W, 2200 m a.s.l), the 
telescope has a 17-m diameter tessellated parabolic mirror, 
supported by a light weight carbon fiber frame. It is equipped 
with a high efficiency 576-pixel 3.5° field-of-view photomul¬ 
tiplier camera. The analogue signals are transported via opti¬ 
cal fibers to the trigger electronics and are read out by a 300 
MSamples/s FADC system. 

At La Palma, HESS J1813-178 culminates at about 47° 
zenith angle (ZA). The large ZA implies a high energy 
threshold of about 400 GeV for MAGIC observations. It 
also provides a large effective collection area (see e.g. 
iKononelko et all ( 1999 ) ). The sky region around the location 
of HESS J1813-178 has a relatively high and non-uniform 
level of light. Within a distance of 1° from HESS J1813-178, 
there are no stars brighter than 8 th magnitude, with the star 
field being brighter in the region SW of the source. 

The MAGIC observations were carried out in the false- 
source tracking (wobble) mode (tFomin et al.ll 19941) . The sky 
directions (Wl, W2) to be tracked are chosen such that in the 
camera the sky field relative to the source position is similar to 
the sky field relative to the mirror source position (anti-source 
position). The source direction is in both cases 0.4° offset 
from the camera center. In Fig. [Tjthese two tracking positions 
are shown by white stars. During wobble mode data taking, 
50% of the data is taken at Wl and 50% at W2, switching 
(wobbling) between the 2 directions every 30 minutes. This 
observation mode allowed a reliable background estimation 
least affected by the large ZA and inhomogeneous star field. 

3. DATA ANALYSIS 

HESS J1813-178 was observed for a total of 25 hours in 
the period June-July 2005 (ZA < 52°). About 15 million 
triggers have been recorded. The calibration of the raw data of 
the MAGIC camera is explained in lGaug et alJ 120051 . Image 
cleaning tail cuts were applied: Pixels are only considered 



Fig. 2.— Distributions of 0 2 values for the source and anti-source, see 
text, for SIZE > 600 ph. el. (corresponding to an energy threshold of about 
1 TeV). 

to be part of the image if their reconstructed charge signal 
is larger than 10 photo electrons (ph.el.) (core pixels) or if 
their charge is larger than 5 ph. el. and they have at least one 
neighboring core pixel. These tail cuts are accordingly scaled 
for the larger outer pixels of the MAGIC camera. The camera 
images are characterized by image parameters (jHilias§T~985l) . 
After the image cleaning and rejection of broken runs about 
10 million events remained for further analysis. These data 
were processed for ' /hadron separation, in a similar way as 
described in iFeganlill9971) . 

In th is analysis, the Random Forest method (see lBock et alJ 
J 200 l for a detailed description) was applied for the 7 /hadron 
separation and the energy estimation. For the training of 
the Random Forest a sample of Monte Carlo (MC) gener¬ 
ated 7 -showers was used to represent the 7 -rays and a ran¬ 
domly chosen sub-set of the measured data was used to rep¬ 
resent the background. The MC 7 -showers were generated 
between 47° and 54° ZA with energies between 10 GeV and 
30 TeV. The spectral index of the generated differential spec¬ 
trum dN^/dE ~ E r was chosen as T = —2.6, in agreement 
with the MAGIC-observed energy spectrum of the Crab neb¬ 
ula dWagner et al.lll200 5). 22 The source-position independent 
image parameters SIZE, WIDTH, LENGTH, CONC IhOdI 
11985 ) and the third moment along the major image a xis, as 
well as the source-position dependent parameter DIST ( Hi 11 as! 
1 9851) . were selected to parameterize the shower images. Af¬ 
ter the training, the Random Forest method allows to calculate 
for every event a parameter, dubbed hadronness, which is a 
measure of the probability that the event belongs to the back¬ 
ground. The 7 -sample is defined by selecting showers with a 
hadronness below a specified value. An independent sample 
of MC 7 -showers was used to determine the efficiency of the 
applied cuts. 

For each event, its o riginal sky position is de termined by us - 
ing the DISP-method dFomin et alll994llLessard et alJ200H) . 
At this stage only source independent image parameters are 
used in the RF training. Figure^ shows the sky map of 7 -ray 
candidate events (background subtracted) from the direction 
of HESS J1813-178. It is smoothed with a two-dimensional 
Gaussian with a standard deviation of 0.1°. To provide a good 
angular resolution a tight hadronness cut and a lower SIZE cut 
of 600 ph. el. has been applied. The SIZE cut corresponds 

22 In order to develop and verify the analysis at high zenith angles, Crab 
data in the interesting ZA range around 50° were taken in lanuary 2005. 
From that sample, we determined the Crab energy spectrum and found it to 
be consistent with other existing measurements (see Fig.^] upper curve). 
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FIG. 3.— Reconstructed VHE 7 -ray spectrum of HESS J1813-178. The 
spectral index is —2.1 zb 0.2 and the integral flux above 400 GeV is about 8 % 
of the Crab nebula (statistical errors only). The dashed line shows th e spec¬ 
trum of the Crab nebula as measured by MAGIC IWagner et al|200^. The 
dot-dashed line shows the results of the HESS collaboration lAharomanet ali 

HSU). 

to an energy threshold of about 1 TeV. The sky map is over- 
layed with contours of 90 cm VLA radio and ASCA X-ray 
data from lBrogan et al . ( 2005 1. The excess is centered at (RA, 
DEC)=( 18 h 13 m 27 s , -17°48'40") and coincides well with the 
position of SNR G12.82-0.02. The systematic pointing un¬ 
certainty is estimated to be 2' (for description of the MAGIC 
telescope drive system see iBretz et aD ( 2003 V) and might in 
future be greatly reduced with the MAGIC starfield monitor 
dRiegel et alJlt2005 ). Apart from the main excess coincident 
with HESS J1813-178 there are no other significant excesses 
present. 

Figure |3 shows the distribution of the squared angular dis¬ 
tance, 6 2 , between the reconstructed shower direction and the 
nominal object position. The observed excess in the direction 
of HESS J1813-178 has a significance of 10.6 standard devi¬ 
ations (according to equation 17 of iLi&Mal fl98l T Within 
errors the source position and the flux level are compatible 
with the measurement of HESS lAharonian et al.l2005blt) . 

For the determination of the energy spectrum, the RF was 
trained including the source dependent image parameter DIST 
with respect to the nominal excess position. A loose cut on the 
hadronness was used. Above the low energy turn-on, the cut 
efficiency reaches about 70% corresponding to an effective 
collection area for 7 -ray showers of about 180000 m 2 . Figure 
0 shows the reconstructed very high energy 7 -ray spectrum 
of HESS J1813-178 after the unfolding with the instrumental 
energy resolution. The measured spectral points are fit by a 
simple power law spectrum, taking the full instrumental en¬ 
ergy re solution into account as described in Mizobuc hi et all 
(2005). The result is given by (x 2 /n.d.f = 5.3/5): 

^ = (3.3 ±0.5) x 10 _ 12 (-E/TeV ) _2 ' 1±o ' 2 

dAataE 

cm _ 2 s” 1 TeV^ 1 . 

The quoted errors (lcr) are purely statistical. The systematic 
error is estimated to be 35% in the flux level determination 
and 0.2 in the spectral index. Within errors the flux is steady 
in the time scales explored within these observations (weeks), 
as well as in the year-long time-span between the MAGIC and 
HESS pointings. 

4. DISCUSSION 

Shortly after the discovery of HESS J1813-178, X-ray 
emission was found in ASCA data coming from the source 



Fig. 4.— Leptonic and hadronic models for the J1813-178 data. Details 
are given in the text. Radio data is from VLA, Bonn, Parkes, and Nobeyama 
observatories (Brogan et al. 2005), X-ray and hard X-ray data are from ASCA 
and INTEGRAL (Ubertini et al. 2005, Brogan et al. 2005). 

AX J 1813- 178 (a.k.a. AGPS 273.4-17.8), iBrogan et all 
d2005l) . also lUbertini et all d2005l) . X-ray emission detected 
by ASCA is predominantly non-thermal, and compatible with 
that expected from a PWN or a SNR shell. X-ray pulsed emis¬ 
sion has not been detected, but the quality and amount of the 
data does not imply strong constraints ferogan et al . 2005 ). 
Statistically, X-ray data are not good enough to unambigu¬ 
ously separate a pure power law from power law + thermal 
contribution either. However, in both cases, none of the two 
component fits yield a significantly different absorbing col¬ 
umn density or photon index when compared with a single 
power law fit dBrogan et al. ) 2005 ). Data analysis from the 
INTEGRAL satellite also showed a soft and luminous source 
at the same location, in the 20-100 keV range (Ubertini et all 
120051) . 

Simultaneously with this X-ray match, HESS J1813-178 
was also found as a non-thermal source in radio data, using 
observations with VLA (90 and 20 cm), Bonn (3 cm) , Parkes 
(6 cm), and Nobeya ma (3 cm) telescopes (Brogan et al 120051: 
iHelfand et al.ll2005i) . These groups discovered a shell-type 
supernova remnant (SNR G12.8-0.0) with a section of the 
shell coinciding with HESS J1813-178. The radio spec¬ 
tral index was found to be -0.48 ± 0.03. There are no 
known radio pulsars detected at the HESS J1813-178 position 
(|Manchester i etd|2005|). 

[Brogan et all (l2005i) conclude that SNR G12.8-0.0 should 
lie at or beyond the distance of W33, ~ 4 kpc. They have 
derived a high column density of Nh ~ 10 23 cm -2 from the 
ASCA data which suggests a significant source of absorption 
in the foreground. 

The multi-wavelength emission coming from the direction 
of HESS J1813-178 is shown in Fig. ^ including the new 
MAGIC data at high energies. We have compared hadronic 
(neutral pion decay) and leptonic (inverse Compton) emis¬ 
sion models with the high energy 7 -ray data, for a review 
see [Torres et all 12003 ). In the case of hadronic models, the 
observed 7 -ray luminosity (2.5 x 10 34 erg s _1 between 0.4 
- 6 TeV, at 4 kpc) implies that the required density of mat¬ 
ter is ~ 6 cm -3 . The 7 -ray production region is presumed 
to be the whole SNR volume. An acceleration efficiency of 
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hadrons of the order of 3% and a supernova power of 10 51 
erg was assumed. Therefore, relativistic hadrons need only 
about 2 Mq of target mass inside the SNR volume (~ 1.5 
pc radius at 4 kpc distance) in order to be consistent with 
the observed luminosity. Alternatively, the target mass for 
the cosmic ray spectrum can be provided by a small cloud 
of a few solar masses located in a region close to the SNR 
shell. For the model shown in Fig. we have assumed 
that the proton distribution is described by dN p /(dVdE ) = 
A p (E/GeV)~ a exp (— E/E max ) GeV - 1 cm -3 where A p is 
a dimensionless normalization constant. We found that a = 
2.1 and E mSLX = 100 TeV is a good fit to the data, and that 
the normalization constant is such that the amount of super¬ 
nova explosion energy converted into relativistic cosmic rays 
needs not to be more than a few percent to agree with MAGIC 
observations. 

For leptonic models, we have assumed a similarly de¬ 
scribed distribution of relativistic electrons dN e /(dVdE) = 
A e (E/ GeV) -ae exp (— E/E max ^) GeV - 1 cm -3 . We found 
that several different inverse Compton models produce rea¬ 
sonable good fits at high energies, e.g. a e = 2.0 and 
-E’max.e = 20 TeV; and a e = 2.1 and 2.2 and B ma x,e = 30 
TeV, all having their energetics at ease with the energy as¬ 
sumed to be released by the supernova explosion. The source 
of target photons for these models was assumed to be the cos¬ 
mic microwave background. The radio spectrum at lower en¬ 
ergies is fitted best by a slope of a e = 2.0. It is computed as 
synchrotron emission of the same electron population. In the 
model shown in Fig. 0] a magnetic filling fraction of about 
20%, a magnetic field of 10 pG and E maxe = 20 TeV and 
a e = 2.0 have been adopted. Thisjnodel is similar to one 
of the models presented bv lBrogan et al.l 12005 ) (blue line in 
their Fig. 3, E maxe = 30, a e = 2.0), without yet having 
a high energy 7 -ray spectrum. The soft X-ray data (ASCA) 
plotted correspond to those shown by Brogan et al. (2005) 
(black data points in their Fig. 3). Our model, as well as Bro¬ 
gan et alls, is in rough agreement with MAGIC data, radio 
and roughly also with the X-ray data, even if the spectral be¬ 
haviour in the ASCA energy range (the slope) is somewhat 
different. That is, we have priviledged the radio data in fitting 
this SED as an example, a better X-ray fitting could be ob¬ 
tained at the expense of a worse radio one. It is worth notic¬ 
ing th at the hard X-ray data of INTEGRAL JUbertini et al.l 
12005 ) cannot be fitted well with the same electron popula¬ 
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tion that is assumed to produce higher and lower energy pho¬ 
tons. Nevertheless, as one can see in Fig. 2 of ( Uhertini et al J 
2005), the ASCA, MAGIC and HESS sources match spatially 
well, while the INTEGRAL/IBIS source is spatially only 
marginally compatible with them. Further multi-frequency 
observations in the x-ray and hard x-ray regime are needed 
to obtain definite conclusions up to what level the leptonic 
models need two electron populations. 

5. CONCLUDING REMARKS 

The detection of HESS J1813-178 using the MAGIC Tele¬ 
scope confirms a new very high-energy 7 -ray source in the 
Galactic Plane. A reasonably large data set was collected from 
observations at large zenith angles to infer the spectrum of this 
source up to energies of about 10 TeV. Between 400 GeV and 
10 TeV the differential energy spectrum can be fitted with a 
power law of slope F = —2.1 ± 0.2. These data can be used 
to cross-calibrate HESS and MAGIC, their independent ob¬ 
servations show satisfactory agreement. 

Multifrequency data in the radio. X-ray, and 7 -ray 
band imply a con nection between HESS J1813-178 and 
SNR G12.82 -0.02 (Helfa nd et al.ll2005t lutertini et al.ll2005t 

Broga n et al.ll2005l) . Generally, hard 7 -ray spectra are ex¬ 

pected from SNRs due to Fermi acceleration of c osmic rays 
iGinzburg & SvrovatskiilHl964l iTorres et al.ll2003h . The hard 
spectrum determined for HESS J1813-178 may be a further 
hint for its association with the SNR G12.82-0.02. 

Present data are not sufficient to discriminate between exist¬ 
ing models for different acceleration mechanisms. Future ob¬ 
servations at lower energies with improved gamma-ray tele¬ 
scopes and/or the GLAST satellite will undoubtedly permit to 
shed more light on the existing leptonic and hadronic models. 
Decisive information concerning hadronic acceleration mech¬ 
anisms is also likely to come from future neutrino telescopes 
like IceCube (lAhrens et al . 200 l . 
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